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Q J25 Recent results of the KASCADE-Grande experiment provided evidence for a mild knee-like struc- 

ture in the all-particle spectrum of cosmic rays at E = p) 16 ' 92±010 eV, which was found to be due 

[/j ■ to a steepening in the flux of heavy primary particles. The spectrum of the combined components 

of light and intermediate masses was found to be compatible with a single power law in the energy 
range from 10 16 ' 3 eV to 10 1S eV. In this paper, we present an update of this analysis by using data 
with increased statistics, originating both from a larger data set including more recent measure- 
ments and by using a larger fiducial area. In addition, optimized selection criteria for enhancing 
light primaries are applied. We find a spectral feature for light elements, namely a hardening at 
E = io 1708±008 eV with a change of the power law index from -3.25 ± 0.05 to -2.79 ± 0.08. 
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High energy cosmic rays exhibit a power law like energy e.g. the EAS-TOP installation [J,|j]. In case of a rigid- 
spectrum with some features, which have been subject to ity dependent steepening in the spectra of the different 
investigations for many years. For energies above 10 eV 50 primaries, the knee in the spectrum of the heavy compo- 
the flux gets too low to be measured directly, therefore nent should appear at an energy of around 10 17 eV jjjlg]. 
one has to extract information about the primary par- According to results from the KASCADE-Grande cxpcr- 
ticles from extensive air showers (EAS), for example, by iment [7| a knee- like feature in the spectrum of heavy 
measuring the secondary particles with ground-based de- particles is visible at E = io 16 - 92±004 eV, as recently 
tector arrays. The most prominent features are the steep- 55 published [8(. It was shown that indeed the bending in 
ening (knee) of the spectrum at about 4 x 10 15 eV and its the all-particle spectrum at E = io 16 - 92±0 - 10 eV is caused 
recovery to a harder slope (ankle) at about 4 x 10 18 eV. by a decrease in the flux of the heavy component. The 
The knee is due to a decrease in the flux of light primaries, combined spectrum of light and intermediate primaries, 
as shown by the KASCADE experiment [l|,l2| and others, however, was therein found to be compatible with a single 

60 power law. 

It is generally argued (see e.g. [3, [K|) that the transi- 
tion from galactic to extragalactic origin of cosmic ra ys i s 
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FIG. 1. (color online). Layout of the KASCADE-Grande 
experiment. The dashed line marks the fiducial area used 
in [g|. The area used in this analysis is indicated by the 
dotted line. The experiment is located at KIT in Karlsruhe. 
The Grande stations (rectangles) had to be arranged to fit 
between the buildings. Therefore the stations are irregularly 
distributed and the missing station at around (—600, —100) 
is the reason for the rectangular cut between —80 m < Y < 
—200 m. The original KASCADE array is located in the upper 
right corner. A comprehensive description of the detector and 
the reconstruction procedures can be found in [7J|. 



tion dominated by light primaries. The present study 
aims to search for experimental evidence for such a spec- 105 
tral feature. We use the same reconstruction and analysis 
procedures as described in [8|, but with increased statis- 
tics and optimized selection criteria for enhancing light 
primaries. 

The number of charged particles, N c h, and the num- 
ber of muons, N^, are the obscrvables used in this anal- 
ysis. Both values are defined as total numbers of parti- 
cles reaching the observation level, where a lateral dis- 
tribution function is fitted to the measured densities at 
the detector stations and integrated. The layout of the 115 
KASCADE-Grande experiment is shown in Fig.[TJ where 
the Grande array records the number of charged parti- 
cles exceeding an energy of 3 MeV by means of 37 scin- 
tillation detector stations, distributed over an area of 
700 x 700 m 2 . The total number of muons is determined 120 
by the measurements of shielded scintillation detectors 
at the KASCADE array, where the threshold energy for 
vertically incident muons is 230 MeV. A detailed descrip- 
tion of the event reconstruction is given in Q . 

A first step to increase statistics is to define a larger 
fiducial area than that used in the first analysis [8[ (see 
Fig. [T]) . Since the probability to misreconstruct the num- 
ber of muons increases with larger distances to the KAS- 
CADE detectors and for lower energies, events with core 
positions farther away from the KASCADE array are 



omitted resulting in the selected fiducial area. By in- 
creasing the fiducial area by 28 %, we also shift the en- 
ergy threshold to higher energies. This is no problem for 
the present analysis, because we are mainly interested 
in the development of the spectrum of light elements at 
energies well above threshold. After including additional 
87 complete days of data taking, the total increase in 
statistics is about 36 % above energy of full efficiency. 

An adequate combination of N c ^ and N^ is used to es- 
timate the primary energy. To classify the events into 
primary mass groups, we use the ratio of N c h to N^ 
in terms of the parameter k (Eq. {]"}). As a function of 
log 10 7V c h, k is centered around for protons and it in- 
creases with the mass of the primary particle to become 
1 for iron. In the same way, the reconstructed energy for 
a given log 10 N^ increases with k from the proton energy 
towards the energy for iron primaries. The k parameter 
is then used in the energy-calibration- function (Eq.Q, 
see also reference fl2l|). 



log 10 (AWAV) - log 10 (iVch/AV)H- 
log 10 (A^ ch /Ag Fo - log 10 (AT cll /Ag H ' 



(1) 



log 10 (iV ch /A''^)H,Fc = C H ,Fc ' lOg 10 (7V ch ) + d H ,Fe, (2) 



log 10 (-E/GeV) = (a H + (a Fo - oh) 



k) • log 10 (AT ch ) 
(b Fc -b H )-k (3) 



flH.Fo, ^H.Fe, ch,Fg and dn.F C are obtained by fitting lin- 
ear functions to the mean log 10 -Etme (coefficients a and 
b) and to the mean \og 10 (N c i 1 /N^) (coefficients c and d) 
of simulated events as a function of the logarithm of 
their reconstructed number of charged particles. This 
has been done separately for 5 zenith angle intervals of 
equal exposure with the upper limits 16.7°, 24°, 29.9°, 
35.1° and 40°, to take the shower attenuation into ac- 
count. The parameters are slightly different from those 
used in |8| , because they are determined using the larger 
fiducial area and taking into account the higher energy 
threshold. For the simulations of the air showers, the 
CORSIKA code [H] was used employing the QGSJET- 
II-2 high energy hadronic interaction model 14 1. Fluka 
(version 2002.4) [15( was used for interactions at low en- 
ergies. The systematic/statistical uncertainties in deter- 
mining the all-particle spectrum as well as the energy 
resolution of the present analysis are shown in Table HI 



TABLE I. Energy resolution and systematic and statistical 
uncertainties for the all-particle spectrum are given for three 
different energies. A more detailed discussion of the uncer- 
tainties can be found in reference lal. 





10 ie6 eV 


10 17 - eV 


10 178 eV 


syst. uncertainty 


11.6% 


11.1% 


24.4 % 


statistical error 


0.5% 


1.2% 


12.0 % 


energy resolution 


22.1 % 


20.5 % 


16.1 % 



1.4 


~ 


^ Data 

. H 






i He 






' CNO 


1. 2 














• Si 




1 
















■ Fe 


1 


• 




i 


< 






1 


* 


* 


4 


4 


« 




OR 


_ — ^ 












y 
















i 












































0.6 
04 


* 


i 


* 


* 


i 


^■T-l 


-i-_T ■ 


£ " 






p 




'"~t- 


0.2 



: 1 


1 

t 




! 




-0.2 


— separation heavy/light between CNO and Si > 

- separation heavy/light between hie and CNO 








-0.4 


_ band of systematic uncertainty 


1,,, 



16.6 16.8 17 17.2 17.4 17.6 17.8 18 

log 10 (E/eV) 









• 


simulated all-particle spectrum - 


- QGSJet 










• 


simulated all-particle spectrum - 


-Epos 










* 


simulated H 


+ He spectrum 


- QGSJet 




10 16 


• 




• 


• 


• • 










- ■ 




• ■ ; 


i 


! ■ 


• 

T 

■ 


• 
* 


• 
■ 


10 1S 


— a 






□ 


m 


□ 








~~ ■ 


rec 


. light spectrum - 


QGSJet 


□ 




5 


□ 




* 


simulated H spectrum - EPOS 








A 




D 


rec 


. light spectrum - 


EPOS 

I , , 






I , , 





16.6 



16.8 



17 



17.2 



17.4 17.6 

log i0 (E/eV) 



FIG. 2. (color online). The mean values of k over the recon- 
structed energy are shown for events simulated with QGSJet- 
II-2, and with zenith angles between 0° —24° , and five different 
primaries. The distributions of k over E are similar in adja- 
cent zenith angle intervals. By combining two angular bins, 
we increase simulation statistics. The increase of systematic 
uncertainties is taken into account. For comparison, the mea- 
sured data (empty squares) is also shown (shifted from the 
bin center to the right for better visibility of the error bars). 
The error bars represent the RMS for the measured data and 
the error of the mean for the simulated data. The dotted 
line was used in [g] to separate the events into a light and a 
heavy component. For an enhanced light selection on cosmic 
rays, the continuous line is taken for the separation, where the 
dashed lines depict the uncertainty of the separation, taking 
into account also the reconstruction uncertainty of k. 



In the following, k is used to separate the all-particle 
energy spectrum into two mass groups: the heavy 
(electron-poor) and the light (electron-rich) component. 
This is done by fitting a separation line to the simu- 
lated mean fc-values over the reconstructed energy and 1! 
comparing the calculated fc-valuc of each event with the 
value given by the separation line. The distributions of 
k for 5 different masses are shown in Fig. [2] Shower- 
to-showcr fluctuations are larger for lighter primaries. 15 
The dotted line was used in [8( to separate the events 
into mass groups. The line is obtained by fitting the 
k Bcp (E) = [fccNo(-E')/2 + k Si (E)/2] distribution, i.e. be- 
tween the CNO and Si simulations. In the present analy- 
sis a lower fc-cut is used to enhance the light component, i< 
which consists of events with assigned fc-values smaller 
than the values shown as solid line (obtained by a fit to 
k scp (E) = [k Hc (E)/2 + k CN o(E)/2}). The dashed lines, 
finally, represent the estimate of a possible error of the 
separation, taking also the reconstruction uncertainty ofi< 
k into account. They are obtained by shifting k scp (E) 
up/down by the statistical and systematic uncertainties 
of k before the fit is performed. A more detailed de- 
scription of the energy determination procedures and the 



FIG. 3. (color online). Shown are simulated and recon- 
structed all-particle energy spectra using QGSJet-II and 
EPOS-1.99 as hadronic interaction models, where for EPOS a 
smaller number of events are available. For the simulations, a 
composition of five elements (H, He, CNO, Si, and Fe) equally 
abundant has been used. The combined proton and helium 
component for QGSJet-II and the pure hydrogen component 
for EPOS are also displayed. In addition, reconstructed spec- 
tra for the light component (separation between He and CNO) 
are shown, where in both cases the QGSJet based reconstruc- 
tion and selection criteria are applied. In case of QGSJet-II 
the light component is reproduced, whereas for EPOS simula- 
tions the reconstructed light component shows an agreement 
with the initial pure proton spectrum. 



separation into mass groups can be found in |8[ and [lfj ■ 
Mainly due to the shower-to-shower fluctuations, it is 
not possible to derive spectra for individual elements us- 
ing k on an event-by-event basis. However, it should be 
noted that in the relevant energy range these fluctuations 
are nearly energy independent. By separating the events 
into light and heavy mass groups according to the present 
procedure, we can assume that these components contain 
at least most of the light/heavy individual elements, re- 
spectively, according to the hadronic interaction model 
employed. 

A simple attempt to cross-check the model dependence 
is displayed in Fig.[3J where we reconstruct events gener- 
ated with EPOS (version 1.99 [l7|) applying the QGSJet 
calibration. The spectra have been simulated with a 
slope index of -2 and weighted to E~ 3 to better repre- 
sent the index of the measured data in this energy range. 
For the simulations, a composition of five elements (H, 
He, CNO, Si, and Fe) with equal abundances has been 
used. The reconstructed light spectra show a significant 
difference in composition, where EPOS generated data 
result in a much lighter composition. This is probably 
caused by the fact that EPOS predicts more muons com- 
pared to QGSJet-II and, therefore, the ratio of A c h to 
N^ is smaller for a given number of charged particles 
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FIG. 4. (color online). The all-particle and electron-rich 
spectra from the analysis [g| in comparison to the results of 
this analysis with higher statistics. In addition to the light 
and heavy spectrum based on the separation between He and 
CNO, the light spectrum based on the separation on He is 
also shown. The error bars show the statistical uncertainties. 



FIG. 5. (color online). The reconstructed energy spectrum 
of the light mass component of cosmic rays. The number 
of events per energy bin is indicated as well as the range of 
systematic uncertainty. The error bars show the statistical 
uncertainties. 



resulting in a larger k- value. Especially helium events 
migrate (by calibrating with QGSJct-II) to the heavy 
mass group. This effect might be slightly compensated 
by the higher reconstructed energy of the events [18[ • Us- 
ing an EPOS calibration, the measured showers appear 
to originate from lighter primaries and of lower energy 
compared to the QGSJet-II calibration. Figure [3] also 
demonstrates that the selection of events according to 
the fc-paramcter does not induce any artificial structures 
in the spectra of light primaries. If the data are well 
described by QGSJct-II, then the spectrum of light pri- 
maries with the separation between He and CNO should 
consist mainly of protons and helium, maybe with some 
additional, less abundant elements between helium and 
carbon. This can be seen in Fig. [3J where the combined 
simulated proton and helium component for QGSJet-II 
is in good agreement with the reconstructed spectrum of 
light elements, which has been obtained by applying the 
QGSJet-II based reconstruction and selection criteria to 
the data simulated using QGSJet-II. Assuming that the 
data simulated with EPOS are closer to real data, then 
the measured spectrum of light particles is an almost 
pure proton spectrum. The simulated proton spectrum 
for EPOS is similar to the reconstructed spectrum of light 
primaries, which has been derived from EPOS generated 
events using again the QGSJct-II based reconstruction 
and selection criteria. According to QGSJet-II, the spec- 
trum of heavy elements for the same separation would 
contain carbon and primaries heavier than that. For 
EPOS it should also contain most of the helium com- 
ponent. 

In Fig. 21 the results of the present analysis are shown. 
To cross-check the results from [8| the all-particle spec-205 
trum and the spectrum of light primaries for the former 



used area and data are compared with the ones obtained 
with higher statistics from the present studies. Both all- 
particle spectra and spectra of light elements based on 
the separation between CNO and Si are in good agree- 
ment. The spectra of light and heavy particles with 
the separation between He and CNO are obtained us- 
ing the separation-line shown in Fig. [5J The spectrum 
of the heavy component, which now contains also the 
medium mass component, exhibits a change of index at 
E = io 16 - 88±003 cV and it therefore agrees inside the 
corresponding uncertainty with the previous result [8| at 
eV. The hardening or ankle-like 
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feature visible in the enriched spectrum of light primaries 
is more prominent compared to the one that includes the 
CNO component. Although statistics gets quite low for 
the spectrum of light elements with the separation on He 
(obtained by a fit to the mean k- values for He in Fig. [2]), 
it is obvious that it cannot be described by one single 
power law only. Formula (QJ [I9[ is used for fitting the 
spectra of the light and heavy components: 



dJ E 

—( E ) = i Q .E^.[i + (—y]^-^ 

lo : normalization factor, 
71/2 : index before/after the bending, 
E'b : energy of the break position, 
e : smoothness of the break. 



(4) 



As shown in Fig.[5j a change of the spectral index from 
71 = -3.25 ± 0.05 to 72 = -2.79 ± 0.08 at an energy of 
±q17.08±o.08 eV is observed for the light component. The 
dashed lines mark the systematic error band for the sep- 
aration between He and CNO obtained by using the se- 
lection shown in Fig. [5J The measured number of events 



above the bending is N meas = 595. Without the bending 
we would expect iV exp = 467 events above this ankle- 
like feature. The Poisson probability to measure at least245 
^Vmeas events above the bending, if -/V exp events are ex- 

pected, is P(N > N meas ) = EZn^J^^'^^) « 
7.23 x 10 -09 . This corresponds to a significance of 5.8 er 
that in this energy range the spectrum of light primaries 
cannot be described by a single power law. If we shift 
the separation criteria in order to obtain an even purer 
proton sample (sep. on He, fig. [4]) the slope difference 
even increases ( 7i = -3.32 ± 0.03 to 72 = -2.59 ± 0.28 
at an energy of io 1716±019 eV). 

So far, possible uncertainties due to the underlying 25 5 
hadronic interaction model have not been taken into ac- 
count. But, as discussed in [20] , applying different meth- 
ods in energy calibration does not result in different con- 
clusions on changes of the spectral slopes. Spectral slopes 
do not significantly change, when corrections for the bin- 260 
to-bin fluctuations (i.e. by unfolding methods correcting 
the reconstruction uncertainties) are applied [16j . There- 
fore, it is not expected that such a correction would result 
in a spectrum of light primaries compatible with a single 
power law. 

One important observation is that the knee in the 
heavy component occurs at a lower energy compared t0265 
the bending in the spectrum of light primaries. This is 
still the case, if the break positions arc shifted by their 
uncertainties. Therefore, the steepening of the heavy 
spectrum and the recovery of the light component is not 
due to a bias in the reconstruction or separation proce-270 
dures. In addition, such a bias should result in a bend- 
ing in the spectrum of light elements at a lower energy 
compared to the spectrum of heavy particles: a bias in 
the separation would be visible in the spectrum of light 
particles first, because its flux is much smaller. A bias275 
in the reconstruction, where the calculated fc-value is too 



low, would result in a reconstruction of heavy particles as 
light particles with a lower energy. It is worth to point 
out that the slope of the heavy mass spectrum above 
the knee-like feature is very similar to the slope of the 
light mass spectrum before the ankle-like feature. This 
could be a hint to a similar source population, as de- 
scribed e.g. in |21| . Correspondingly, the slope index of 
72 = —2.7 might be an indication of an injection spec- 
trum of a new (extragalactic) source population of high 
energy cosmic rays [9h11|. 

In summary, after redefining what is considered as 
a electron-poor/electron-rich event and increasing the 
statistics, a significant hardening, i.e. an ankle-like fea- 
ture in the cosmic ray spectrum of light primaries, is 
observed at E = io 1708±008 eV. The slope index of 
the underlying power law changes at this energy from 
-3.25 ± 0.05 to -2.79 ± 0.08, which might be an indi- 
cation that the transition from galactic to extragalactic 
origin of cosmic rays starts already in this energy range. 
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